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ABSTRACT 

We present Spitzer infrared observations of the starless core L429. The IR images of this core show 
an absorption feature, caused by the dense core material, at wavelengths < 70 /im. The core has a 
steep density profile, and reaches Ay > 35 mag near the center. We show that L429 is either collapsing 
or in a near-collapse state. 

Subject headings: ISM: globules - ISM: individual (LDN 429) - infrared: ISM - (ISM:) dust, extinction 



1. INTRODUCTION 

The collaps e of a dense cloud core into a star hap- 
pens rapidly (|Havashilll966l ) and therefore is rarely ob- 
served; furthermore, most dense cores appea r not to be 
collap sing but to be close to equilibrium (e.g-. lLada et al.l 
2008). If they are not sufficiently supported by gas 
pressure and turbulence, then a plausible level of mag- 
netic field strength would suffice to support them (e.g ., 
iKandori et alll2005t IStutz et al.ll2007t lAlves et al.ll2008j ) . 
Because the process of collapse is so rarely observed, 
there is considerable debate about what triggers it and 
how it proceeds to form a star. Additional examples of 
dense cores in extreme states that indicate either collapse 
or a condition just prior to collapse are critical to con- 
strain theories about the first stages of star formation. 

Mm-wave line observations are the traditional means 
to search for coll apsing cloud cores, s uch as in the pio- 

^ (119861). The subsequent 
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discussio n (e.g. , iMenten et al.1 



Woottenl 119891 : IMundv et al 



19981) illustrates the ambiguities in the interpretation of 
velocity data in star forming dense cores. Alternative ob- 
servational approaches are necessary to supplement such 
line measurements and provide a better understanding 
of collapsing and nearly-collapsing dense cores. 

In this paper we characterize the structure of the 
extremely dense core L429, using the infrared shadow 
techn ique we have demonstr ated on two other glob- 
ules (jStutz et al.l I2007L l2008f l This core is already 
known to have an extremely high deuterium fraction, 
believed to be an indication of its being very close to 
collapsing into a protostar, and a low internal temper- 
ature of only Tki n ~ 8 K, reducing pressure suppor t 
to low levels (IBacmann et alJl20"o"l ICrapsi et aLll2005[ K 
iBacmann et al.l ( 20001 ) detected a 7 //m shadow from 
L429 using ISOCAM. Here, we report that this shadow 
is so opaque that it is detectable at 70 /im. The presence 
of significant absorption at this long wavelength demon- 
strates that L429 is indeed extremely dense. Our analysis 
of the shadow also demonstrates that it is very compact 
with a steep density gradient. We show how these char- 
acteristics put L429 in the collapsing or near-collapsing 
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starless core category. 

2. OBSERVATIONS AND PROCESSING 

L429 was observed by Spitzer using MIPS, program 
ID 30384, P.I. Tyler L. Bourke, and with the IRAC in- 
strument, program ID 20386, P.I. Philip C. Myers. The 
MIPS observations utilized scan map mode in all three 
b ands. The data red uctions were carried out as described 
in IStutz et all (j2007f ). The MIPS images are - 15f5 x 54' 
in size, with exposure times of ~ 150 s at 24 /im, ^60 s 
at 70 /im, and ~ 15 s at 160 /im. For the 70 /im and 
24 /im data respectively the beam size is ~ 18" and ~6", 
the nominal default pixel scale is 9'.'8 and 2'.'5, and the 
pixel scale in our images is 4" and 1'.'24. The IRAC re- 
ductio ns were carried out as described in lHarvev et al.l 
(2006]). The IRAC images are ~6' x 5f5 in size, with ex- 
posure times of ~ 800 s. The Spitzer data are shown 
in Figure Q] We also use the 2MASS All- Sky Point 
Source Catalog (PSC) (|Skrutskie et all 12006) photome- 
try for two well-detected sources within 40" of the center 
of L429Reddening estimates for these sources, assuming 
the iRiekeT fe Lcbofsky (1985) extinction levels, are plot- 
ted in Figure O 

3. THE DISTANCE 

Th e commonly accepted distance to L429 is 200 pc 
fe-g.. iLee et~aHll999b IBacmann et~aLll2000t ICrapsi et all 
l2005f ). based on t he likely association of thi s cloud with 
the Aquila Rift. iDame fc Thaddeud (|1985h find a dis- 
tance to the Aquila Rift based on the CO line Vlsr — 
8 km s _1 . They summarize the results of various studies 
and conclude that the distance to the Aqu ila Rift is con- 
sistent with 200 ± 100 pc. More recently, IStraizvs et al.l 
(2003) derive a di stance to the nea r side of the Rift of 
225 ± 55 pc, while Kawamu ra et al.l (|2001l ) derive a dis- 
tance to an Aquila cloud complex centered at ~ I = 48° 
and r^b = —6° of 230 ± 30 pc, both consistent with pre- 
vious studies. The depth of the absorption feature (see 
Fig-ID implies that L429 must lie on the near side of, or in 
front of, the Aquila Rift because substantial amounts of 
foreground material would be warmed by the interstellar 
radiation field and would dilute the shadow. We assume 
a distance of 200 pc, consistent with previous distance 
determinations. If the distance to L429 is larger the de- 
rived mass will increase, making the core less stable and 
more likely to be collapsing. 

4. SHADOW PROFILES 




Fig. 1 — MIPS and IRAC 5' X 5' images of L429 centered at RA = 18 h 17 m 05.4 s , Dec = -08°13'32.3" (J2000; I = 21.68°, b = 3.80°). 
The circle (radius = 40") indicates the location of the shadow produced by the dense core material. The images are displayed on a log 
scale; north is up and east is to the left. 



L429 is seen in absorption at A < 70 /im; at 160 /im 
there is large-scale emission but no isolated local peak is 
observed (see Fig. [1]). In addition, no compact emission 
source is detected at any Spitzer wavelength within the 
shadow, and so no source of luminosity > 10~ 3 L is 
present (jHarvev et all 120061 : [Dunham et alj I2008T ). We 
derive extinction profiles from the observed shadows at 
70, 25, 8.0, and 5.8 /n m using techniques simil ar to those 
discussed in detail in lStutz et all (|2007l 120081 ). 

4.1. 70 fim profile 

The first step in deriving the 70 /im profile is to sub- 
tract from the image a value that represents zodiacal 
emission and other foreground sources. For L429 the 
shadow is the darkest part of the image. For this reason 
the image zero level cannot be determined directly from 
the image itself. To set an upper limit on the optical 
depth, we zero the image using the lowest pixel value in 
a 44" x 44" box centered on the shadow; in this case the 
zero value is ~ 5.2 MJy sr _1 below the mean pixel flux 
level immediately outside the shadow. Ideally, one would 
compute the minimum value based on more pixels. How- 
ever, the 70 /xm pixel size is a non-negligible fraction of 
the size of the core of absorption feature; taking more 
pixels would dilute the depth of the shadow. We use 
IRAS data to set the lower limit on the optical depth 
profile: we measure the flux level outside the shadow to 
be ~ 13.5 MJy sr" 1 above the foreground component of 
emission, interpolated to 70 /xm from the 60 and 100 /im 
IRAS images and estimated far away from regions with 
significant emission (about 4.5° away). For this value, 
the minimum shadow value is 8.3 MJy sr _1 above the 
image zero level. This zero level is a very conservative es- 
timate because the shadow region is surrounded by high 
levels of emission which are very likely to raise the fore- 
ground contribution above the level we use here. This 
lower limit also accounts for th e MIPS 70 /im calibration 
error, 5% ([Gordon et al.ll2007D . and the IRAS extended 
source calibration erro r, which can be as large as 30% 
([Wheelock et al.lll99l . 

After zeroing the image we mask out all pixels with 
possible emission above the sky level. The mask thresh- 
old is set to the mean value plus 0.5c measured for pixels 
in a 4' x 4' region 4' north of the shadow. Next we mea- 
sure the mean pixel values, termed / here, in nested, and 
concentric, circular regions, with radii starting at 7'/5 and 
increasing in 5" steps, out to ~ 120", centered at RA — 
18 ft 17 m 05.4 s , Dec = -08°13'32.3" (J2000), without ac- 



counting for fractions of pixels. We ignore all pixels in 
the southern half of the shadow because the shadow has a 
tail extending to the south and the assumption of spher- 
ical symmetry breaks down (see Figure [lj . Finally, to 
determine the optical depth at 70 /im, T70 = —ln(f/fo), 
we average profile values at radii > 27" 5 to measure /o, 
the intrinsic, unobscured local flux value. The outer ra- 
dius of the core at 70 /im, equal to 27.5", is set by the 
observed flattening of the measured flux profile at that 
radius, indicating that there is litt l e to n o absorption at 
larger radii. Following IStutz et~a l. (2007), we derive the 
70 /im Ay profile shown in Figure [2] with open trian- 
gles: downward-pointing triangles represent the upper 
limit and the upward-pointing triangles represent the 
lower limit. The innermost region has an optical depth 
of T70 = 1.8 to 0.46 (upper and lower limits). The total 
mass (gas plus dust) implied by the 70 /im profile is given 

by 



M 



70 



5.1 M, 







/ D \ 2 1.58 x 10 2 cm 2 gnr 1 

V200 pC/ Kabs,70 



(1) 



where we have assumed a gas to du st mass ratio of 100 
and the lOssenkopf &; Hennin g (1994) model dust opacity 
for coagulated dust grains of 1.58 x 10 2 cm 2 gm" 1 . The 
value of 5.1 Mq is the upper limit on the mass; the lower 
limit mass is 1.6 Mq. The enclosed mass profile is shown 
in Fig. [3] As noted earlier, our limits on r are very 
conservative. We performed an error analysis to measure 
the effect of a reasonable perturbation in the assumed 
zero level and the calculated value for /o and find that 
the errors are absorbed into the limits by a wide margin. 
For a dust opacity of Av a bs. 7o = 6.35 x 10 1 cm 2 gm -1 , from 
the Ry = 5.5 model from IDraind (|2003allbT l. we calculate 
a 70 /im mass of M 70 = 12.6 - 4.0 Mq. 

4.2. 24, 8.0, and 5.8 /im profiles 

The method used to derive the shorter wavelength pro- 
files is essentially the same as described above. First, we 
zero the images using the darkest pixel in the shadow. 
Second, the mask threshold is determined from the pixel 
value distribution in the sky region. We measure the 
mean pixel value in circular regions starting with a ra- 
dius of 2"5 and increasing in steps of 2"5, out to a radius 
~ 120". The value of /o is determined by averaging pixel 
values between radii of 40" and 45". These radii were 
chosen based on a flattening of the 24 /im profile in this 
region. At radii between 45" and ~ 55" the profile con- 
tinues to rise, likely due to contamination from other 
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Fig. 2. — Extinction profiles for L429. The downward- and 
upward-pointing triangles (connected with a dotted line) represent 
upper and lower bounds on the 70 fim extinction profile, respec- 
tively; the innermost 70 (im region has an optical depth of T70 = 1.8 
to 0.5, for the upper and lower bounds, respectively. The dashed 
lines represent the Bonnor-Ebert best-fit models for the 70 /an ex- 
tinction profile derived assuming a distance of 200 pc; both models 
require high temperatures that are in contradiction with the ob- 
servations (see § 5). The 24 fim (squares), 8.0 fim (diamonds), 
and 5.8 fim (plus signs) profiles are scaled to match the outer two 
radial regions of the 70 fim profile; for these three profiles we only 
show points with r < 1.0. The crosses (x) are 2MASS extinction 
measurements for two source s with reliable detections. We assume 
Osscnkopf & Hcnning ( 1994) dust opacities for the extinction pro- 
files (see § 4). 



nearby sources. Using larger values of /o would artifi- 
cially increase r. This choice of /o is further justified 
by the 8.0 and 5.8 /im profiles: they show an increase 
in the mean pixel value well past 40" but a decreasing 
value approaching the inner estimate at radii larger than 
^80", and by a radius of ~ 100" the profiles match the 
inner values. This increase in the profile values at in- 
termediate radii may be caused by limb brightening of 
the core as it is irradiated by the interstellar radiation 
field. If the emission raises the observed shadow flux then 
it will cause an underestimate of r at all wavelengths; 
however, if it predominantly increases the unobscured 
flux levels outside the shadow r will be over-estimated. 
We do not attempt to correct for this effect. We cal- 
culate optical depth limits of T24 = 2.1,r 8 . = 3.2 and 
T5.8 = 2.2 through the center of the profiles. The im- 
plied total masses at these wavelengths is ~ 1.2 for 
all three profiles, assuming Osscnko pf fc Hennind (|1994f ) 
model dust opacities. The discrepancy between this mass 
and the 70 /im mass (between 1.6 and 5.1 M Q ) most 
likely results from the high optical depth in the core. If 
the core is dense enough to be observed in absorption at 
70 /im then it is certainly opaque at shorter wavelengths; 
the Ay profile and mass estimates at shorter wavelengths 
will be lower limits. In the CB190 case the core was ob- 
served in absorption at 24 fim but not at 70 /im, and the 
maximum measured extinction value was Ay ~ 30 mag 
(T24 = 1.4), much smaller than the values measured here 
jstutz et al J 120071) . Therefore we reject all 24, 8.0, and 
5.8 /im r > 1.0 points and we scale the profiles to match 
the 70 /im extinction values at the outer two 70 /im radii. 
In Figure Owe plot the scaled profiles; the shapes of the 
profiles match well. 

5. DISCUSSION 



° 2 



D = 200 pc 

V 



□ □ □ □ 



VA 70 fim 
□ 24 fim 



10 



20 

radius (arcsec) 



30 



40 



Fig. 3. — Enclosed mass vs. radius for L429 calculated assuming a 
distance of 200 pc. The downward- and upward-pointing triangles 
represent masses calculated using the upper and lower bounds on 
the 70 fim extinction profiles, respectively. The 24 fim masses 
(squares) are scaled to match the 70 fim profile; we only show 
24 fim points with r < 1.0. We assume Osscnkopf & Hcnning 
(1994) dust opacities (see §4). 



Bonnor-Ebert models are often used to character- 
ize starless cores; these models generally fit the ob- 
served density profi les with reasonable fi delity (e.g., 
i Kandori et ail 120051 : iShirlev et ail 120051 : IStutz et all 
l2007f h Bonnor-Ebert models describe a self-gravitating 
pressure confined isothermal sphere in hydrostatic equi- 
librium. The density profiles are described by three 
parameters: temperature, central density p c (or equiv- 
alently, £ max ), and the outer radius of the cloud. A 
model will be unstable to collapse if £ max > 6.5. We fit 
said models to the extinction profiles discussed above. 
We use a 3D grid of models in temperature, radius, 



and £ max ; see IStutz et ail (|2007l ) for a more detailed 
discussion of the models and the fitting method. One 
strong constraint in interpreting the best-fit models is 
that L429 is dark at A < 70 /im, and therefore the 
temperature cannot be greater than ~ 12 K, in agree- 
ment with other core temperature measu r ements (e.g., 
i Lemme et al.l 119961: iBacmann et al.l l2002t iHotzel et ail 
l2002tlKirk et al.ll2007l ). The best-fit models for the up- 
per and lower bounds on the 70 /im profile are shown 
in Figure [Has grey dashed lines. The 70 /im profile has 
a best fit temperature of T = 56 — 130 K, for the up- 
per and lower bounds, respectively; these temperatures 
are in clear contradiction with the data. The best-fit 
outer radius is 27", for both the upper and the lower 
limit, and the best-fit £ max = 4.5 and 3.0, for the upper 
and lower limit respectively. As the assumed distance is 
increased (in this case to 270 pc) the best-fit models be- 
come even more implausible, with temperatures always 
> 70 K. At a distance of 100 pc the best-fit temperatures 
are T = 26 — 90 K, still in contradiction with the data. 
These high temperatures indicate that L429 requires ad- 
ditional support (turbulence and magnetic fields) beyond 
t hermal to pre v ent co llapse. 

ICrapsi et al.1 (|2005l ) observe N 2 D+ and N 2 H+ line 
widths ~0.3 — 0.4 km s _1 and calculate excitation tem- 
peratures T cx ~4.5 — 5 K. The particle of mean molecular 
mass has velocity FWHM ~0.4 km s _1 for temperature 
8 K, indicating that the turbulent energy contribution to 
the support in L429 is only comparable to the thermal 
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contribution (jHotzel et al.ll2002f) . 

Magnetic fields are a possible source of support in star- 
less cores. Inferred magnetic field s trengths in cores 
range from, e.g., < 15 ^G to 160 /iG (|Bergin fc Tafallal 
120071 ). For 9 = 30", a distance of 200 pc, and B = 
160 uG, the mass that can be supported is Mb = 1 M Q 
(Stahl eTfe Pallal I2005T). Furth ermore, from equation 1 
and the lstahler fc Pallal (|2005f ) analysis, the condition of 
stability requires a minimum magnetic field strength of: 



B = 861 [iG 



1.58 x 10 2 cm 2 gm" 1 



^abs,70 
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30' 



where 9 is the angular radius of the core. The mean 
column density of the 70 fim profile is (N(R 2 )) - 1.1 X 
10 23 cm" 2 ; giv en this column, the correlation presented 
by iBasul (|2004f ) implies a line-of-sight magnetic field of 
Bios = 250 /iG; we estimate the total field as B tot = 
v3 x Bi os — 430 /iG, although this relation has large 
scatter. 

Although the field required to support L429 is larger 
than expected, our understanding of the limits on mag- 
netic field strengths does not allow for a strong con- 
straint. In addition to requiring a field exceeding the 
likely range of values, supporting the globule entirely 
with magnetic pressure has other consequences. One 
might expect asymmetries to develop, given the asym- 
metric pressure provi ded by magnetis m. If such struc- 
tures can be avoided (jShu et al.lll987f ). another issue is 
that the field should leak out of the glo bule by ambipo- 
lar di ffusion, on timescales ~3 x 10 6 yr (|Stahler &: Pallal 
2005). Thus, either the field is already too small to pre- 
vent the collapse of the globule, or the collapse will only 
be delayed for a time. 

If infall is occurring in L429, th e observed line width is 
very narrow. ICrapsi et alJ (120051) detec t a double-peaked 
N 2 H+ (1 - 0) line and lLee et alJ (|1999ft detect a double- 
peaked CS (2— 1) line, both measurements wit h roughly 
equal brightness in the blue and red peaks. iLee et all 
(1999) state that because other thin tracers show sin- 
gle Gaussian components, the observed double peak is 
likely due to self absorption at high optical depth; these 



authors also note that L429 is one of two sour ces out of 
163 sp ectra that show small wing components. ILee et ahl 
( 2004) classify L429 as a possible infall candidate based 
on their observation of CS (3 — 2), which has a dou- 
ble peaked profile with more emission in the blue peak. 
iSohn et all (|2007l ) also measure a double peaked profile 
in HCN (1 — 0); they note that this source is one of two 
(out of 85 cores surveyed) that has some line strength 
anomalies in the th ree hyperfine HCN lines. Recently, 
ICaselli et all (|2008fl have found that L429 has a large 
amount of ortho-H2D + (li > o-li,i) emission, which they 
find to be correlated with increased central concentra- 
tion, CO depletion, and deuteration, providing further 
evidence that L429 is very evolved. 

6. CONCLUSIONS 

We measure the density profile of L429 from the 70 fim 
absorption feature observed with the Spitzer Telescope. 
We conclude that: 

1. The strong absorption even at 70 /Ltm indicates that 
L429 is exceptionally dense. 

2. At the same time, the absorbing cloud is very com- 
pact and has a very steep density gradient (as seen most 
clearly at 24 fim). 

3. The temperature can be no more than 12 K. 

4. As a result of these extreme characteristics, thermal 
pressure fails, by a wide margin, to support the globule; 
turbulent support is also inadequate by a wide margin. 

5. The magnetic field required to make up the rest of 
the required pressure for support would be surprisingly 
large; even if present, such a large field will leak out of 
the globule through ambipolar diffusion. 

6. Therefore, L429 is either already undergoing collapse 
or is approaching an unstable, near-collapse state. 
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